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Calculations of the potential energy surface for tracer Ga and In adatoms above three GaAs 
(111) A surface reconstructions are presented in order to understand the growth conditions required 
to form axial heterostructures in GaAs/InGaAs nano-pillars. In all calculations the Ga adatom has 
a stronger bond energy to the surface than the In adatom. The diffusion barriers for Ga adatoms 
are 140meV larger than for In adatoms on the Ga vacancy surface, but they are comparable on the 
As trimer surface. Also the binding energy for an In adatom is closer to that of a Ga adatom on 
the As trimer surface. We conclude that the As trimer surface is preferable for adsorption of In and 
thus for selective formation of hetero-interfaces on (111) facets. This work helps explain the recent 
successful formation of axial GaAs/InGaAs hetero-interfaces in catalyst free nano-pillars. 



INTRODUCTION 

Semiconductor nanowires(NWs) and nanopillars(NPs) 
are exciting materials for probing mesoscopic physics 
and as building blocks for future high performance opto- 
electronic devices on Si 1 3 . NP synthesis by catalyst-free 
selective area metal-organic chemical vapor deposition 
(SA-MOCVD) is a growth technique for forming large 
arrays of uniform NPs in lithographically defined loca- 
tions with the inclusion of optical alignment marks for 
device integration 4 . 

The absence of a metal particle to catalyze growth 
means that atoms adsorb directly onto the crystal sur- 
faces from the vapor, and the resulting crystal shape is 
controlled in part by minimization of the total surface 
energy 5 . GaAs nanopillars grow in the [111] direction, 
and have hexagonal symmetry with side facets composed 
of the {011} family of planes. Atoms from the vapor 
adsorb on all facets of the NP and then diffuse to the 
(111) surface at the tip where they incorporate. The po- 
lar (111) surface has a higher surface energy than the 
stoichiometric {011} planes, making the observed crystal 
shape energetically favorable. 

Heterostructure formation is a necessary capability to 
master in catalyst-free NP synthesis in order to create ef- 
ficient optical device^. Core-shell hetero-structures have 
been studied in a variety of material systems, but ax- 
ial hetero-structure formation has been elusive in this 
growth mode. When a new atomic species is introduced, 
the surface energetics must promote incorporation of the 
new species on the top (111) surface while simultaneously 
suppressing incorporation on the side walls. Despite this 
challenge, axial InGaAs segments of varying composition 
and thickness were recently demonstrated in GaAs cata- 
lyst free NPs grown by SA-MOCVE0. High V/III ratios 
(V/III ~ 50) were required to promote incorporation 



of In in the axial direction with negligible shell growth. 
At the lower V/III ratios (V/III ~ 10) typically used 
for GaAs NP homoepitaxy, indium is not selective to the 
(111) surface, and instead nucleates on the side- walls, de- 
forming the crystal facets. Fig.[TJa) shows scanning elec- 
tron micrograph (SEM) of NPs formed by SA-MOCVD 
with axial InGaAs inserts at high V/III ratio, and the 
vertical side- walls and hexagonal symmetry are evident. 
Fig. [TJb) shows a dark field scanning transmission elec- 
tron micrograph (STEM) of the same pillars revealing 
the axial InGaAs segment. In contrast, Fig. [TJc) shows 
pillars terminated with InGaAs sections at low V/III ra- 
tios, and having deformed crystal facets due to indium 
nucleation on the side- walls. This tendency for indium 
to bond to all available crystal surfaces has also been 
reported in Ref [8]. 



To investigate possible reasons for the observed dif- 
ferences in behavior between In and Ga during nanopil- 
lar epitaxy, we present a theoretical investigation of the 
potential energy surface (PES) for Ga and In tracer 
adatoms situated above three common surface recon- 
structions of GaAs(lll)A. The technique of calculating a 
PES has been applied by numerous researchers as a tool 
for studying diffusion, adsorption and desorption and for 
understanding epitaxy on crystal surfaces^ffiJ] A simi- 
lar study of In and Ga tracer diffusion on GaAs {011} 
is necessary for a more complete understanding of NP 
epitaxy, and will be presented in a future publication. 
Computational methods are discussed first, followed by 
a description of the calculations and their results. We 
conclude with a discussion and interpretation of the re- 
sults. 
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FIG. 1. (a) SEM of GaAs nanopillars containing axial InGaAs 
inserts grown at high V/III ratio, (b) Dark field STEM of sin- 
gle InGaAs insert, (c) SEM of GaAs nanopillars terminated 
with InGaAs at low V/III ratio. 



COMPUTATIONAL METHODS 

To calculate the potential energy surface (PES) of a 
Ga or In adatom above a GaAs(lll)A surface recon- 
struction, we begin by computing the equilibrium surface 
geometry of the three reconstructions depicted in Fig. |2j 
From left to right, the surfaces are the Ga vacancy sur- 
face, the As trimer surface and the As adatom surface. 
All three surfaces possess a 2x2 unit cell indicated by a 
shaded parallelogram. Slabs 9 mono-layers thick are it- 
eratively relaxed, keeping the bottom three mono-layers 
fixed, until residual atomic forces are < 0.02eV/A. 

The total energy of the surface with an additional Ga 
or In adatom is then computed using a 4x4 super cell. 
The entire system, slab and adatom, is allowed to relax, 
but the adatom coordinates are fixed perpendicular to 
the [111] direction (the adatom is fixed in the x-y plane 
and allowed to relax in z). All three surfaces possess 
3-fold rotational symmetry, and each rotationally sym- 
metric slice posses a mirror symmetry such that only 6-8 
points are sampled in a triangle above the 2x2 unit cell. 
The calculated energies are then reflected, rotated twice 
through 120° and mapped to a rectilinear grid using a 
cubic interpolation to generate a PES for the adatom of 
interest. The energy zero-point is chosen to be the total 
energy of the relaxed reconstructed surface plus the total 
energy of an isolated atom of In or Ga. 

Calculations were performed within the framework of 
density-functional theory (DFT) as implemented in the 
software package FHI-AIMS 12 , which uses numeric atom 
centered orbit als for its basis set. The Perdew-Burke- 
Ernzerhof (PBE) parameterization of the generalized gra- 
dient approximation is used for the exchange correlation 
functional. Approximately 16 layers of vacuum and 64 
equivalent k-points in the lxl unit cell are specified. Con- 
vergence of the energy difference between the maximum 
and minimum on the PES is confirmed for the k-points, 
slab thickness, vacuum layers and super-cell size for the 
Ga vacancy reconstruction. In addition, total energy dif- 
ferences were tested for the FHI-AIMS built in settings 
"light" and "tight" 12 . The "light" setting, having fewer 
basis functions and a smaller numerical integration mesh, 
is found to increase the speed of the calculation while gen- 



TABLE I. Diffusion barrier Ed, minimum potential energy 
A±, secondary minimum potential energy A2 and transition 
point T of In and Ga adatoms above three common surface 
reconstructions of GaAs (111) A. All values are in electron 
volts (eV). 
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Ga Vacancy 


Ga 


1.06 


-2.87 


-2.21 


-1.81 


-1.73 
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-2.65 
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As Trimer 


Ga 
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-6.83 
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0.26 


-6.99 


-6.88 


-6.73 




As Adatom 


Ga 


0.20 


-6.88 




-6.68 






In 


0.13 


-6.71 




-6.58 





erating results that differ from "tight" by only a few meV. 
Calculations are therefore carried out using the "light" 
setting. 



RESULTS 

The potential energy surfaces for indium and gallium 
adatoms above each surface reconstruction are presented 
in this section. The binding energies at adsorption sites, 
Aj, and transition points, T and X", for In and Ga 
above each surface are collected in Table H The domi- 
nant diffusion energy barriers, calculated as the differ- 
ence Ed =T — Ai, are also tabulated. The Ga vacancy 
and the As trimer surfaces are the primary surfaces of in- 
terest because they are energetically favorable in a vapor 
consisting of mixed As and Ga atoms^H The Ga vacancy 
surface has the lower surface energy at low As chemical 
potentials and the As trimer surface as the lower surface 
energy at high As chemical potentials. The As adatom 
reconstruction always has the highest relative surface en- 
ergy, and is presented here for completeness, even though 
this surface does not exist with high probability. 

The PES for Ga and In adatoms above the Ga vacancy 
surface reconstruction are shown in Fig. |3| This surface, 
characterized by a missing Ga atom, is thought to be 
the stable reconstruction under most conditions except 
in extremely arsenic rich environments 1 ^. Comparing In 
and Ga adatoms above the Ga vacancy surface, the PES 
are qualitatively similar, with a deep minimum at the 
vacancy site Ai, and a secondary minimum at the site 
A2, above third layer As atoms. The transition points, 
T and T' , are saddle points of the PES that are crossed 
when hopping between adsorption sites. 

Diffusion can occur on the Ga vacancy surface by two 
possible pathways. Either the adatom hops directly be- 
tween A\ sites over the transition point X", or it crosses 
over the point T into the secondary site A2, and then 
rapidly hops back into an adjacent A\ site. The energy 
barrier T' — A\ is 80 meV larger than the transition over 
X. At typical growth temperatures of ~ 1000 dif- 
fusion between A\ sites by way of A2 is fast enough to 
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FIG. 2. Surface reconstructions of the GaAs(lll)A surface. From left to right the Ga vacancy surface, the As Trimer surface 
and the As adatom surface. Arsenic atoms are depicted as light gray spheres, and gallium atoms are depicted as dark gray 
spheres. Top and side views are rendered with two or three layers of atoms The atomic diameters are drawn larger for atoms 
closer to the surface. The 2x2 unit cell is identified by a shaded parallelogram. While all surfaces possess a 2x2 unit cell, the 
PES calculations were performed using a 4x4 super-cell to suppress interaction effects between the adatoms. 



dominate the diffusion path. The diffusion barrier, Ed, 
reported in table [T] is the barrier to hop from A\ to A^. 

Ga atoms are less mobile than In on this surface with a 
diffusion barrier 140 meV higher than In regardless of the 
path taken. The binding energy of a Ga adatom at A\ is 
220 meV larger than for In, suggesting that Ga adatoms 
will be adsorbed preferentially over In adatoms. This 
calculation agrees with the observation that In floats to 
the surface when forming the NP hetero-interface^. 

The PES for a Ga adatom above the Ga vacancy re- 
construction was previously calculated by Taguchi et. al. 
in Ref. [10]; however, our results are significantly differ- 
ent. In that work, contrary to expectations, they found 
the potential energy minimum was not in the lattice site 
vacated by the Ga atom, but at adjacent interstitial lo- 
cations with diffusion energy barriers of ~ 0.4eV. Our 
calculations, in contrast, show a deep potential minimum 
at the vacant lattice site with diffusion barriers ~ l.OeV. 
We are unable to explain the discrepancies between the 
two calculations, however, the authors of Ref [10] ac- 
knowledge that the vacant Ga site should be more stable 
according to calculations based on the interatomic po- 
tential. In light of the conflicting results, we carefully 
checked our energy calculations and algorithms for gen- 
erating the PES, which exploit the surface symmetry, and 
are unable to find errors in our methods. 

The As trimer PES for In and Ga adatoms are pre- 
sented in Fig. [4] The As trimer surface is the stable 
reconstruction appearing in arsenic rich environments. 
The PES for both In and Ga adatoms have potential en- 
ergy minimum A\ at the center of the As trimer, and a 
diffusion barrier height of 260-270 meV. The PES for an 
In adatom also has a secondary minimum, A2 , above one 
of the second layer As atoms that can potentially slow 
the diffusion for In. The difference in binding energy be- 
tween Ga and In is only 110 meV for the As trimer sur- 
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FIG. 3. Gallium vacancy potential energy (PES) surface for 
Ga (left) and In (right) adatoms. The top atomic layer of the 
reconstruction is drawn as an overlay to assist in visualizing 
the adsorption sites. The primary potential minimum occurs 
at Ai with a secondary minimum at A2. At typical NP growth 
temperatures, diffusion is dominated by hops between A\ and 
A2 over the saddle points labeled T. 



face, compared to 220 meV for the Ga vacancy surface. 
Indium adatoms will have a higher probability of adsorp- 
tion on this surface compared to the Ga vacancy surface 
because of the equivalent diffusion coefficients and more 
competitive binding energy. 

The PES for the In and Ga adatoms above the As 
adatom surface are shown in Fig. [5] The stable adsorp- 
tion sites, Ai, are above the top layer As atoms, but the 
entire region surrounding the As adatom is energetically 
very flat with minor variations of tens of meV. The re- 
gion surrounding the As adatom is therefore an adsorp- 
tion site. Diffusion barriers are lowest for this surface 
reconstruction compared to the other surface reconstruc- 
tions with a barrier of 130 meV for In and 200 meV for 
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FIG. 4. Potential energy surface for Ga (left) and In (right) 
adatoms above the As Trimer reconstruction of GaAs(lll). 
The top atomic layers of the reconstruction are drawn as an 
overlay to assist in visualizing the adsorption sites. 
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FIG. 5. Potential energy surface for Ga (left) and In(right) 
adatoms above the As Adatom reconstruction of GaAs(lll). 



Ga. The absolute binding energy is 170meV larger for 
Ga than for In on the As adatom surface reconstruction. 



DISCUSSION 

Calculations were performed to provide a physical ex- 
planation for why high V/III ratios are required for for- 
mation of GaAs/InGaAs axial heterostructures in (111) 
oriented catalyst free NPs. We believe that the high As 
chemical potential results in surface reconstructions on 
the NP that suppress In incorporation on the {011} NP 
side-walls while simultaneously promoting In incorpora- 
tion on the (111) NP tip. 



The calculations reported in this work support the hy- 
pothesis that the As trimer surface, stable at high As 
chemical potential, is desirable for higher rates of in- 
dium incorporation on the (111) surface for two reasons. 
First, the difference in binding energy between Ga and 
In adatoms in the A\ adsorption site is reduced from 
220meV on the Ga vacancy surface to HOmeV on the As 
trimer surface. This reduction means that In adatoms 
compete more effectively with Ga and have a higher prob- 
ability incorporating into the crystal in the presence of 
the As trimer surface. Second the diffusion barriers, Ejj 
are comparable for both Ga and In adatoms on the As 
trimer surface, yet the diffusion coefficient of In is roughly 
two orders of magnitude larger on the Ga vacancy sur- 
face at typical growth temperatures of ~ 1000i\~. In the 
presence of a Ga vacancy surface, In adatoms will diffuse 
more quickly than Ga and desorb from the small (111) 
surface at the tip of the pillar. The resulting chemical 
environment of adsorbates at the pillar tip will be richer 
in Ga than in the surrounding vapor. In contrast, the 
comparable diffusion barriers of both Ga and In on the 
As trimer surface will result in a concentration of ad- 
sorbates representative of the concentration in the sur- 
rounding vapor. The two reasons cited explain why In 
adatoms incorporate more efficiently on the (111) surface 
at high As chemical potential. 

In summary, calculations of the PES for tracer In and 
Ga adatoms above three stable surface reconstructions 
of GaAs(lll)A were performed. Gallium has a larger 
binding energy to all surfaces, yet the binding energy of 
In is competitive with Ga on the As trimer surface. In 
addition, the diffusion barriers for In and Ga are iden- 
tical on the As trimer surface, which forms at high As 
chemical potentials. On the Ga vacancy surface, stable 
at low As chemical potentials, indium is calculated to 
rapidly diffuse, lowering the residence times at adsorption 
sites, reducing the opportunity for incorporation into the 
crystal. These results suggest that formation of an As 
trimer surface at high V/III ratios promotes formation 
of axial GaAs/InGaAs hetero-interfaces during nanopil- 
lar growth. 
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